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ABSTRACT 

We perforin a joint analysis of X-ray and Sunyaev Zerdovich (SZ) effect data using an analytic model 
that describes the gas properties of galaxy clusters. The joint analysis allows the measurement of the 
cluster gas mass fraction profile and Hubble constant independent of cosmological parameters. Weak 
cosmological priors are used to calculate the overdensity radius within which the gas mass fractions 
are reported. Such an analysis can provide direct constraints on the evolution of the cluster gas mass 
fraction with redshift. We validate the model and the joint analysis on high signal-to-noise data from 
the Chandra X-ray Observatory and the Sunyaev-Zel'dovich Array for two clusters, Abell 2631 and 
Abell 2204. 

Subject headings: X-rays: galaxies: clusters-galaxies: individual (Abell 2631, Abell 2204) 



1. INTRODUCTION 

Galaxy clusters trace the growth of structure in the 
universe. Their abundance and evolution is critically sen- 
sitive to underlying cosmological parameters such as Om, 
fTg, and the dark energy equation of state parameter w. 
Recent work has focused on using galaxy clusters to con- 
strain cosmology, including dark energy constraints from 
X-ray and joint X-ray/SZ measurements of the gas mass 
fraction (Allen et al. 2008; Ettori et al. 2009; LaRoque 
et al. 2006), cosmological parameter constraints from the 
growth of structure via X-ray (Mantz et al. 2008, 2010; 
Vikhlinin et al. 2009) and SZ cluster surveys (Vander- 
hnde et al. 2010; Marriage et al. 2011; Sehgal et al. 2010; 
Muchovej et al. 2011; Williamson et al. 2011; Benson 
et al. 2011). 

In this paper, we present a method for the joint anal- 
ysis of X-ray and SZ cluster observations using a self- 
consistent analytic model for the physical properties of 
the intra-cluster medium (Bulbul et al. 2010). The model 
provides analytic expressions for the radial density, tcm- 
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perature, and pressure profiles, and is therefore simul- 
taneously applicable to both X-ray and SZ observables. 
The joint analysis allows measurement of the cluster gas 
mass fraction without the need to impose external pri- 
ors on cosmological parameters such as the Hubble ex- 
pansion rate H{z). Such an analysis applied to a sam- 
ple of clusters can directly probe the evolution of clus- 
ter gas mass fractions with redshift. We demonstrate 
the method using high signal-to-noise data from Chan- 
dra and Sunyaev-Zel'dovich Array (SZA) observations of 
two clusters, Abell 2631 and AbeU 2204. 

The method developed in this paper combines all avail- 
able Chandra X-ray data (both imaging and spectro- 
scopic) with SZ observations, using these to determine 
the angular diameter distance and cluster mass with min- 
imal cosmological assumptions. We chose the Bulbul 
et al. (2010) model for this analysis since it describes the 
three thermodynamical cluster properties - density, tem- 
perature, and pressure - with a consistent set of param- 
eters that are both readily interpreted. The temperature 
profile linking density and pressure is both easily calcu- 
lated and observationally motivated. This work differs 
from Mroczkowski et al. (2009) who used SZ and X-ray 
imaging data only with a simplified, core-cut form of the 
Vikhlinin et al. (2006) model to describe the X-ray den- 
sity and the Nagai et al. (2007) parameterization of the 
SZ pressure, and for which the inferred temperature pro- 
file did not reduce to a compact, accessible expression. 

This paper is structured as follows: Section 2 describes 
the data reduction and analysis. Section 3 the modeling 
of X-ray and SZ data, and Section 4 the joint analysis 
of the X-ray and SZ data. We present and discuss our 
conclusions in Section 5. 

2. DATA REDUCTION 

We selected a non-cool core cluster (Abell 2631) and 
a cool-core cluster (Abell 2204) with high quality X-ray 
and SZ observations to demonstrate the method of anal- 
ysis. 
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Fig. 1. — SZ contours overlaid on X-ray false color images of Aboil 2631 {left) and Abell 2204 {right ). The Chandra X-ray surface 
brightness data are from the energy range 0.7-7.0 keV. The color bars reflect the number of counts detected by Chandra in the 0.7-7 keV 
band, with a pixel size of 1.97 arcsec. The SZ data are from the SZA, and the contour levels are (-1-2,-2,-4,-6,-8,...) times the rms noise 
(see Table 1). The FWHM of the synthesized beams for these SZ observations are shown in the lower left corner of each image. 



2.1. Chandra Imaging and Spectroscopy 

The Chandra X-ray data are in the form of event files, 
which we use to generate both images and spectra. Addi- 
tional blank-sky composite event files are used for back- 
ground subtraction. The event files are reduced using 
CIAO 4.3.1 and CALDB 4.3, following the reduction pro- 
cedure described in Bulbul et al. (2010). Details of each 
cluster observation can be found in Table 1. 

The X-ray images in the 0.7-7 keV band are used to 
measure the X-ray surface brightness profile of the clus- 
ter. To subtract the background from the surface bright- 
ness, we rescale the blank-sky image to match the clus- 
ter surface brightness in a peripheral region that is free 
of cluster signal. The peripheral regions we chose are 
at a distance > 500 arcsec from the cluster center for 
Abell 2631 (corresponding to 2.7 Mpc in the standard 
flat ACDM cosmology) and > 550 arcsec (1.7 Mpc) for 
Abell 2204. 

Spectra are extracted in annular regions centered at 
the peak of the X-ray emission. These regions cover an 
area out to the radius where the surface brightness pro- 
file reaches the background, which is near rsoo (the radius 
within which the average density is 500 times the criti- 
cal density at the cluster redshift) for these observations. 
From the blank-sky data, background spectra are also 
extracted and processed. We then rescale the blank-sky 
spectra to match the count rate of the cluster spectra 
in the 9.5-12 keV band. In this band, Chandra has no 
effective area for the detection of photons, and the de- 
tected counts originate from a particle background that 
is time variable. Hickox & Markevitch (2006) showed 
that while the flux within the 2-7 keV and 9.5-12 keV 
energy bands can vary with time, the ratio of the two 
bands remains constant. Subtracting the blank-sky data 
rescaled by the higher-energy band therefore accurately 
removes the background from the lower-energy band. 

After rescaling the blank-sky spectra and removing the 
background from our cluster data, residuals may still be 
present in the soft 0.7-2 keV energy band. These soft 



X-ray residuals may be due to Galactic and extragalac- 
tic emission, and may vary as function of position (e.g., 
Snowden et al. 1997) and time (e.g., Takei et al. 2008). 
For each cluster observation, we use a peripheral region 
that is free of cluster emission — the same region used 
to rescale the background images — to determine whether 
soft residuals are present after the blank-sky background 
has been subtracted. We detect the presence of soft resid- 
uals in both clusters. The residual spectra arc flt us- 
ing a phenomenological model that includes a power law 
and a plasma emission model, and this model is rescaled 
by area and included in the spectral fit for each annu- 
lus (e.g., Snowden et al. 1998; Nevalainen et al. 2005; 
Maughan et al. 2008). 

2.2. Interferometric Observations with the SZA 

The two clusters were observed with the Sunyaev- 
Zel'dovich Array (SZA), an eight-element interferometer 
designed for the detection and imaging of the SZ effect. 
Each antenna in the array is 3.5 m in diameter and has 
a primary beam FWHM of 10.7' at the center frequency 
of the observing band (31 GHz). For these observations, 
six antennas were closely packed together to provide sen- 
sitivity to arcminute-scale SZ signals, and the remain- 
ing two antennas were placed further out to constrain 
the flux contributions from unresolved radio sources, as 
described in Muchovej et al. (2007). The unflagged on- 
source time for Abell 2631 was 16.1 hours and 19.6 hours 
for Abell 2204. The details of the observations are given 
in Table 1, and the radio sources detected in each field 
are listed in Table 2. In the analysis of the cluster SZ ef- 
fect described below, the parameters of the SZ decrement 
and the radio sources are fit simultaneously. 

The SZA data are reduced using a set of routines writ- 
ten in MATLAB"'^'* that constitute a complete pipeline 
for flagging, calibrating, and reducing visibility data. 
The reduction pipeline, described in Muchovej et al. 
(2007), converts the data to physical units and cor- 

http:/ /www. mathworks.com/products/matlab 
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rccts for instrumental phase and amplitude variations. 
Data arc flagged for corruption due to bad weather, 
sources of radio interference, and other instrumental ef- 
fects that could impact their quality. The pipeline out- 
puts calibrated unflagged visibilties, i.e., components of 
the Fourier transform of the sky brightness multiplied by 
the primary beam response, along with their correspond- 
ing statistical weights and positions in the Fourier (u — v) 
plane. 

3. MODELING THE X-RAY AND SZ DATA 

The X-ray observable is spatially resolved spec- 
troscopy, and temperature and metallicity of the intr- 
acluster plasma are measured using the X-ray spectro- 
scopic data. The X-ray surface brightness is defined as 



by Bulbul et al. (2010): 



5x = 



1 



47r(l + z)3 



/ 



nlAee{Te,A)d£ 



(1) 



where £ is the line of sight through the cluster, rie is the 
electron density, Tg is the electron temperature, A is the 
metallicity, and Aee{Te, A) is the X-ray cooling function 
(in units of counts cm^ s~^) as a function of electron 
temperature and metallicity. The density, temperature, 
and metallicity can vary along the line of sight. 

The observable from the SZ data is the amplitude of 
the spectral distortion of the Cosmic Microwave Back- 
ground (CMB) in the direction of the cluster. This dis- 
tortion is due to inverse Compton scattering of CMB 
photons off electrons in the intracluster medium (ICM), 
and results in a decrement in the CMB brightness tem- 
perature at frequencies < 218 GHz. The magnitude of 
the decrement is proportional to the electron pressure 
integrated along the line of sight (Sunyaev & Zel'dovich 
1972): 



(2) 



where Tcmb is the temperature of the CMB, f{x,Te) con- 
tains the frequency dependence of the SZ temperature 
signature using the relativistic corrections provided by 
Itoh et al. (1998) and Nozawa et al. (2006), ctt is the 
Thomson cross section, nif. is the electron mass, and c is 
the speed of light. 

Assuming spherical symmetry, the line of sight inte- 
gration element d£ relates to the angular element (in ra- 
dians) as d£ = Df^d9, where Da is the angular diameter 
distance. Prom Equations 1 and 2, we find 
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The combination of X-ray imaging spectroscopy and SZ 
observations can be used to simultaneously measure the 
distribution of the electron density, the electron temper- 
ature, and the angular diameter distance (e.g., Hughes 
& Birkinshaw 1998; Grego et al. 2000; Reese et al. 2002; 
Grainge et al. 2002; Saunders et al. 2003; Bonamente 
et al. 2006). For further discussion on the SZ effect and 
its use for cosmology, see reviews by Birkinshaw (1999) 
and Carlstrom et al. (2002). 

We describe the density and temperature profiles of the 
hot plasma in galaxy clusters using the model proposed 
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rieo is the normalization of the pressure profile, n is the 
polytropic index, Tg is the scale radius, To is the normal- 
ization factor for the scaling of the temperature profile, 
7 is the slope of the cooling function, and a is the cool- 
ing parameter which ranges from to 1. One attractive 
feature of these models is that they provide a simple an- 
alytic form for the electron pressure: 



Pe(r) =PeO(/)(r,r«,;3) 



n+l 



(9) 



where Pgo = rieokTa is the pressure normalization. The 
free parameters of the model, which are used to jointly 
fit the X-ray and SZ data, are Heo, Tq, rg, rcooi, ct, l3,j,n 
and the distance Da- 

Parameter estimation is done using a Monte Carlo 
Markov chain (MCMC) method described in Bonamente 
et al. (2004). Correlation among model parameters is a 
feature of most analytic models, including the one used in 
this paper. Parameter correlations can result in low ac- 
ceptance rates and thus slow convergence of the Markov 
chains (e.g., Gilks et al. 1996). In our implementation, 
steps of the MCMC are proposed along the set of direc- 
tions which diagonalize the covariance of the posterior 
distribuition, determined via singular-value decomposi- 
tion of an initial test chain, resulting in efficient explo- 
ration of the parameter space (see Appendix). 

3.1. X-ray Data Analysis 

The annular bins in the tcmpcratiire profile (Figure 2) 
were chosen by starting with an initial 10" bin and then 
increasing each bin by 50% of the width of the previous 
bin to give roughly the same counts per bin. Follow- 
ing the analysis of the systematic uncertainties for the 
X-ray data described in Bulbul et al. (2010), we adopt 
a 1% systematic uncertainty on the count rate of each 
bin of the surface brightness profile and a 10% system- 
atic uncertainty on the temperature of each spectral re- 
gion as discussed in Section 4.5.1. Figure 2 shows the 
temperature and surface brightness profiles, along with 
the best-fit models, for the fit to the X-ray data only 
of Abell 2631 and AbeU 2204. Since our model has the 
same parameters for both density and temperature, the 
surface brightness profile carries a larger weight in the 
fit. The model fits are acceptable for both clusters to 
within the plotted errors, which include the systematic 
uncertainties associated with the surface brightness and 
temperature discussed in Section 4.5. For Abell 2631, 
there is insuflacient signal to constrain n and ^ simulta- 
neously. In this case, we fix /3 = 2, which is equivalent 
to assuming that the total mass follows a Navarro et al. 
(1997) profile at large radii (see Section 4.3). 
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Fig. 2. — X-ray surface brightness and temperature profiles for Abell 2631 {left column) and Abell 2204 (right column). Top panels: 
Surface brightness profiles where the black points are derived from the X-ray images, the red line shows the best fit model, and the green 
line is the background level determined from the blank sky observations. The residuals show the fractional difference between the model 
and the data. A 1% systematic uncertainty has been added in quadrature to the datapoints. We plotted the surface brightness profiles 
beyond the fitted region to show the agreement between the cluster emission and background. Bottom panels: Temperature profiles where 
the red line shows the best fit model and the dark, medium, and light grey regions show the 68%, 95.4%, and 99.7% confidence levels 
obtained from the model fits. A 10% systematic uncertainty has been added in quadrature to the temperature bins. 



3.2. SZ Data Analysis 

After removal of corapact radio sources in the cluster 
field, the visibilities V,/(u, v) measured by the SZA can be 
related to the Fourier domain equivalent of the integrated 
Compton-y parameter Y(u,v) (see, e.g., Mroczkowski 
et al. 2009). This is defined as 



Y{u,v) 



g{x)Io ' 



(10) 



where g{x) corrects for the frequency dependence of the 
SZ flux, and Iq = 2{kBTc,j,uf I {hcf is the primary CMB 
intensity normalization. Figure 3 shows Y(u, v) along 
with the best-fit model, for the fit to the SZ data only of 
the two clusters. 



4. JOINT ANALYSIS OF X-RAY AND SZ DATA 

We first perform a consistency check of the determina- 
tion of the pressure profiles from the X-ray and SZ data. 
We then focus on determinations of the angular diameter 
distance and the radial profile of the gas mass fraction 
using the consistent parameterization of density, temper- 
ature and pressure provided by the ICM model for the 
joint analysis of the X-ray and SZ observables. 



4.1. Consistency of X-ray and SZ measurements of the 
electron pressure profiles 

X-ray and SZ observations provide independent mea- 
surements of the radial distribution of the electron pres- 
sure. The X-ray observables are the electron temperature 
and surface brightness, the latter depends on the square 
of the electron density according to Equation 1; the SZ 
observable, on the other hand, is directly proportional 
to the electron pressure integrated along the line of sight 
(Equation 2). 

The two observables can be affected by different 
sources of systematic uncertainties. For example, the 
presence of non-thermal X-ray emission (e.g.. Million & 
Allen 2009; Bonamente et al. 2005; Sarazin & Lieu 1998) 
could result in the increase of the X-ray emission above 
the level of the thermal gas, and radio emission from 
cluster halos (e.g., Brunetti et al. 2007) may partially 
fill the SZ decrement. Another source of systematic un- 
certainty is the assumption of spherical symmetry in the 
analysis (e.g., Sulkanen 1999), which would result in a 
different measurement of the pressure from X-ray and 
SZ observations. A discussion of sources of systematic 
uncertainty in the analysis of X-ray and SZ observations 
is presented in Section 4.5. A comparison of the pressure 
profiles from SZ and X-ray observations is therefore use- 
ful to determine the presence of sources of emission that 
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Fig. 3.— SZ visibility profiles for Abell 2631 (top) and Abell 2204 
(bottom) plotted as a function ofu—v radius (\/«^ +v^). The plots 
show the real components of the measured Y{u, v) along with the 
best fit model. 

can cause differences between the two measurements. 

We perform a joint fit to tlie X-ray and SZ data us- 
ing Equation 9, with the normalization of the SZ pressure 
model independent of the X-ray density and temperature 
normalizations. The common parameters in the models 
{ts, P, and n) are linked between the two datasets, thus 
requiring the X-ray and SZ pressure profiles to have the 
same shape. In this analysis, we adopt the angular di- 
ameter distance appropriate for the cluster redshift in a 
ACDM model with h = 0.73, = 0.27 and VLa = 0.73. 

The pressure inferred from the X-ray and SZ measure- 
ments are within 20% of one another for both clusters, 
consistent with the statistical and systematic effects (Ta- 
bles 3 and 6); results for a larger sample of clusters can 
be found in Bonamente et al. (2011). Note that the mea- 
surement of the ratio of X-ray pressure to the SZ pressure 
depends on the choice of the Hubble parameter, since the 
pressure normalizations are degenerate with the value of 
Da assumed in the analysis. 

4.2. Direct measurement of the angular diameter 

distance 

We also perform a joint X-ray and SZ analysis which 

enables us to place direct constraints on the angular di- 
ameter distance without using priors on the cosmological 
parameters (e.g., Hughes & Birkinshaw 1998; Birkinshaw 
1999; Reese et al. 2002; Bonamente et al. 2006). For this 
analysis, we link the shape parameters and the pressure 
normalizations between the X-ray and SZ data, and al- 
low Da to vary. For Abell 2631 and Abell 2204 we mea- 
sure Da = 798.9 i^gfi Mpc and Da = 575.3 ±|f :| Mpc, 



both values are consistent with those calculated using a 
standard a ACDM cosmology at the l-cr level. The mea- 
surement of the Da for Abell 2204 is also in agreement 
with that of Bonamente et al. (2006). The measurement 
of the angular diameter distance for a given cluster is 
affected by a number of systematic effects (Bonamente 
et al. 2006), and the agreement of Da with the ACDM 
value is expected for a large sample but not necessarily 
for individual clusters, as it is for the two clusters in this 
paper. 

4.3. Radial profiles of the gas mass fraction independent 
of cosmological parameters 

By using direct measurement of Da as described in 
Section 4.2, we can also obtain radial profiles of the gas 
mass, total mass, and the gas mass fraction without the 
need to use priors on the cosmological parameters. The 
gas mass Mgas is computed by integrating the gas density 
profile within the volume. 



gas 



Aw Hemp J ne{r)r^ dr = Aw HempD ^ j ne{0)O^d9, 

(11) 

where fx is the mean molecular weight (calculated assum- 
ing metal abundances of 0.3 solar, Anders & Grevesse 
1989), mp is the proton mass, and dr = Dp^d9. The total 
mass Mtot is computed assuming hydrostatic equilibrium 
between the gravitational mass and the thermal pressure 
of the gas: 



Mtot(r) 



Awpirl 



+ 



1/(1 - /3) - r/r. 



(/3-2) V/3-1 (l + r/r, 



Tcoolir), 

(12) 

where = (ToA:(n + l)(/3 - I)) / [AwG^imprl). 

Figure 4 shows the radial profiles of the gas mass, 
total mass, and gas mass fraction for Abell 2631 and 
Abell 2204. The uncertainties reflect the fact that Da 
is also measured directly from the data, and that no as- 
sumption about the value of the cosmological parameters 
i?o, or Oa was made. 

4.4. Measurement of the gas mass fraction at an 

overdensity radius 

In cosmological applications (e.g., via the distribution 
of /gas with redshift, Allen et al. 2008) /gas is typically 
measured within an overdensity radius. The radius ta 
is defined as the radius within which the average matter 
density of the cluster is A times the critical density of 
the universe at the cluster's redshift: 



3 ^ Mtot (r a) 
"^-fAp,(.)' 



(13) 



where pc{z) = {2,HqE'^{z))/{&wG) is the critical density 
of the universe, Ha is the Hubble constant, and E^{z)= 
fluCi- + zf -h Oa -I- r2k(l + z)^ in the ACDM model. 

The joint X-ray and SZ analysis provides cosmology- 
independent constraints on Da and on the radial pro- 
file of /gas(^) (Sections 4.2 and 4.3). The radius ta 
and therefore all quantities calculated out to this ra- 
dius retain a cosmological dependence through the fac- 
tor Pc{z) appearing in Equation 13. In the following we 
describe a method to marginalize the measurement of 
/gas (''a) over the cosmological parameters, which results 
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Fig. 4. — Top Panels; left Abell 2631, right Abell 2204 '■ Gas mass profiles determined from tiie joint analysis of Chandra/SZA observations 
for Abell 2631 and Abell 2204. Middle Panels: Total mass profiles. Bottom Panels: Gas mass fraction profiles. The dashed lines are 
the 68% confidence level at each radius. Grey areas show the measurements at radii r25oo and rsoo, obtained by marginalization over the 
cosmological parameters. Chandra calibration systematics are included in the measurements. 



in a weak cosmology dependence of our joint measure- 
ments of /gas(rA). 

In what follows, wc adopt the standard flat Friedman- 
Robertson- Walker (FRW) cosmological model, parame- 
terized by Ho and Om, with = 1— fijvi, which is known 
to provide a good fit to other cosmic distance data (Allen 
et al. 2008; Frccdman ct al. 2009; Pcrcival ct al. 2010). 
Within this model, a given value corresponds to a 
curve in the {Hq, Qm) plane described by 

fj ^ ^ r .... 

° Da{i + z)J„ [i]M((f + C)^-i)-i]^/2' ^ ' 

where spatial flatness is assumed. For each step in the 
Markov chain of each cluster, wc use the corresponding 
value of Da to generate a consistent pair of cosmological 
parameters by drawing a random value of from a uni- 
form distribution on [0,1], and calculating the associated 
value of Hq from Equation 14. We then use these param- 
eters when calculating ta and the other fit parameters 



of the chain. This approach has the advantage of using 
minimal prior information on the cosmological parame- 
ters, with free to range between and 1, and Da 
measured directly from the data. This method properly 
accounts for the covariance between ta, /gas (''a), Da 
and the other parameters of the fit. For the flat FRW 
model and the low rcdshifts in question, this procedure 
corresponds, to very good approximation, to the use of 
uniform priors on Hq and flu- We note, however, that 
this correspondence does not necessarily hold for higher 
redshifts {z ~1) or for other cosmological models. The 
measurements of masses and /gas at r25oo and r^ao are 
indicated as grey areas in Figure 4 and are listed in Table 
4. 

For comparison, we also provide the measurements of 
masses and /gas using the X-ray data with a fixed Da 

calculated from the standard flat ACDM cosmology, and 
priors on flM from the WMAP measurements of Ko- 
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matsu et al. (2011). The results arc reported in Table 5. 

4.5. Sources of Systematic Uncertainty 

The sources of systematic uncertainty on the gas mass 
fraction and the integrated X-ray and SZ pressure are 
listed in Table 6. The individual errors are added in 
quadrature to determine the total systematic uncertainty 
on fgas and the integrated pressure. Systematics from 
Chandra instrument calibration and the X-ray back- 
ground are included in the fitting of the data, and the 
value of masses and fgas m Tables 4 and 5 account for 
these systematics. Since this joint analysis method can 
be applied to larger samples of clusters, we also indicate 
whether the impact of each source of systematic uncer- 
tainty is reduced with sample size. 

4.5.1. Instrument Calibration 

Chandra^s ACIS effective area has a spatially depen- 
dent non- uniformity at the level of ±1% and therefore 
we add a ±1% uncertainty to the surface brightness 
data. We also adopt a ±10% uncertainty on the tempera- 
ture measurements to account for uncertainty in the low- 
energy calibration of the effective area (see, e.g., Bulbul 
et al. 2010). These uncertainty estimates for the Chan- 
dra calibration are folded into the mass measurements 
reported in Tables 3, 4 and 5. 

Frequent observations of Mars are used to calibrate the 
SZA absolute flux scale; we employ the Rudy (1987) flux 
model which has an estimated absolute calibration un- 
certainty of ±5%. The stability of the instrumental gain 
is ±3%, as determined from repeated calibrator measure- 
ments in SZA survey fields (Muchovej et al. 2011). The 
absolute calibration and instrumental gain yield a global 
±6% uncertainty on the SZA calibration. We rescaled 
the SZA data by the ±6% uncertainty on the SZA cali- 
bration and compared the measurements to the original 
analysis. We found a ±6% systematic uncertainty on the 
pressure and ±10% systematic uncertainty on the gas 
mass fraction at r25oo and rsoo- The uncertainty associ- 
ated with instrument calibration does not average down 
with sample size. 

4.5.2. Kinetic SZ Effect 

Reese et al. (2002) reports that for a cluster with a 
temperature of 8.0 keV and with a typical velocity along 
the line of sight of 300 km/s (Watkins 1997; Colberg et al. 
2000) the kinetic SZ effect would be ±4% of the thermal 
SZ for 30 GHz observations. Accordingly we use a ±4% 
uncertainty due to the kinetic SZ effect on the gas mass 
fraction and the SZ pressure profiles for measurements at 
''2500 and rsoo- This source of uncertainty averages down 
by a factor of the square root of the size of the sample. 

4.5.3. Radio Source Contamination 

Undetected radio sources not accounted in the model- 
ing could lead to a biased measurement of the SZ decre- 
ment. We use the Faint Images of the Radio Sky at 
Twenty-centimeters (FIRST) database as a reference for 
locating compact radio sources within 10' of the cluster 
center. Most radio sources that will affect the 30 GHz 
data (rms noise 0.25 mjy) will have counterparts in 
the FIRST survey (rms noise of 0.15 mJy at 1.4 GHz); 
inverted spectrum sources may not have counterparts at 



1.4 GHz and will affect our measurement of the SZ decre- 
ment, but they comprise a small fraction of the source 
population (Muchovej et al. 2010, Fig. 3). 

We determine the effect of undetected point sources 
by placing a radio source model at each FIRST source, 
fixing the position and marginalizing over the flux. In 
the pressure model (Equation 7), we fixed the parameters 
Ts, n, /3, and Da and let P^o be free. We compare the 
pressure profiles with the original analysis (see Table 2) 
and find a < 1% difference in the pressure of each cluster. 
Therefore we conservatively apply a 1% uncertainty on 
the pressure and 2% uncertainty on the gas mass fraction 
at r25oo and rsoo- 

4.5.4. Asphericity 

Although we assume a spherical model in our analysis, 
most clusters do not appear to be circular in shape in X- 
ray or radio observations. LaRoque et al. (2006) reports 
a 10 — 20% uncertainty in the measurement of the gas 
mass fraction due to asphericity; therefore we use ±20% 
as a conservative estimate for measurements at r25oo and 
rsoo- This uncertainty also averages down by a factor of 
the square root of the sample size, as shown in Sulkanen 
(1999), provided the selection of the sample is unbiased 
with respect to cluster shape. 

4.5.5. Hydrostatic Equilibrium Assumption 

The assumption of hydrostatic equilibrium at large 
radii results in an underestimate of the total mass. This 
is due to the presence of non-thermal pressure which can 
bias hydrostatic equilibrium measurements of the total 
mass (Evrard et al. 1990). According to Lau et al. (2009) 
the total mass of a relaxed cluster such as Abell 2204 will 
be biased by —6% at r25oo and —8% at rsoo, and for un- 
relaxed systems such as Abell 2631 by —9% at r25oo and 
—11% at rsoo- Therefore we adopt a systematic uncer- 
tainty of —9% and —11% in our error analysis for mea- 
surements at r25oo and rsoo- 

4.5.6. X-ray Background 

The X-ray background is determined from the ACIS 
blank sky composite event file. We normalize the 
blank sky background level for each observation using 
an emission-free region on the AGIS detector. We ad- 
justed the background normalization factor by a factor 
of ±2fT and propagated this through the analysis, and 
found that this produces a ±2% uncertainty on the back- 
ground count rate. This uncertainty affects the surface 
brightness and temperature measurements resulting in 
a ±2% and ±9% uncertainty on the gas mass fraction 
measurements at r2soo and rsooj and a ±2% uncertainty 
on the X-ray pressure profiles at both radii. This uncer- 
tainty averages down by a factor of the square root of 
the sample size. 

4.5.7. Systematics associated with the use of the Bulbul 
et al. (2010) models 

The Bulbul et al. (2010) model assimies a polytropic 
relationship between the ICM density and temperature 
at large radii. To estimate uncertainties associated with 
the polytropic assumption, we compare our X-ray masses 
for the two clusters in our sample with those calculated 
using the Vikhlinin et al. (2006) model, which provides 
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an independent parameterization of the thermodynamic 
quantities. From this comparison, we find that the gas 
mass fraction measurements varies by < 10% at all radii 
between the Bulbul et al. (2010) and Vikhlinin et al. 
(2006) models. We estimate the uncertainty on the cor- 
responding pressure profiles by comparing the integrated 
pressures between the two models, and the comparison 
results in an uncertainty of ±3% at ^2500 and ±5% at 
'"500- We consider these uncertainties as a rough esti- 
mate of the systematics associated with the Bulbul et al. 
(2010) model. 

4.5.8. Helium Sedimentation 

The effect of helium sedimentation may be an addi- 
tional source of systematic uncertainty. In our mea- 
surements we assume that the hydrogen to helium ratio 
is uniform throughout the cluster. However, theoreti- 
cal studies (Fabian & Pringle 1977; Rephaeli 1978) sug- 
gest helium sedimentation effects may affect cluster mass 
measurements. Peng & Nagai (2009) finds that the bias 
in gas mass fraction from the presence of helium sedi- 
mentation is less than 10% at ~ r25oo and negligible at 
rsoo- Accordingly, we estimate a systematic uncertainty 
of 10% at r25oo and <5% at rsoo for the measurement 
of the gas mass fraction. Bulbul et al. (2011) applied 
the Peng & Nagai (2009) helium sedimentation simula- 
tion model to a sample of clusters and demonstrated the 
effects on the gas mass and total mass. The integrated 
pressure is proportional to the gas mass, and we use the 
values from Bulbul et al. (2011) to determine upper lim- 
its to the systematic uncertainty in the measurement of 
pressure of —4% at r25oo and —2% at r5oo- 

5. CONCLUSIONS 

We demonstrate the use of the Bulbul et al. (2010) 
cluster model for simultaneous fitting of X-ray data and 
Sunyaev-Zel'dovich effect data. The model employs a 
compact parameterization that relates the three primary 
thermodynamic quantities by the ideal gas law at all 
radii. We consider X-ray data from Chandra and 30-GHz 
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weakly dependent on f^M- 

We discuss possible sources of systematic errors in the 
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TABLE 1 
Cluster Observations 



Cluster 


z 




RA 


Dec 


Chandra Observations 




SZA Observations 








(10^»cm-2) 


(,J2000) 


(J2000) 


ObsID 


Time 


On-sre Time 


FWHM 


P.A. 


rms Noise"* 














(ks) 


(hrs) 


(arcsec) 


(deg) 


(mJy) 


A2631 


0.27 


3.55 


23:37:40.1 


+00:16:33 


3248/11728 


25.0 


16.1 


152 X 117 


17.2 


0.4 


A2204 


0.15 


5.67 


16:32:47.2 


+05:34:32 


7940 


72.9 


19.6 


157 X 115 


-7.7 


0.4 



''FWHM (full-width at half maximum of the synthesized beam), P.A. (position angle of the synthesized beam), and rms noise are for short 
baselines (< 2feA). 



TABLE 2 
Radio Sources in Cluster Fields 



Cluster 


Pointing 


Center 




30 GH 


z Source 








rms Noise'' 




RA 


Dec 


src 




A5" 


Flux 


FWHM 


P.A. 




(J2000) 


(J2000) 




(arcsec) 


(arcsec) 


(mJy) 


(arcsec) 


(deg) 


(mJy) 


Abell 2631 


23:37:38.8 


+00:16:06.5 


1 


21.3 


36.5 


3.7 


26.5 X 16.6 


42.7 


0.25 








2 


205.0 


-130.0 


0.5 








Abell 2204 


16:32:46.88 


+05:34:32.4 


1 


0.4 


1.2 


7.0 


21.1 X 18.4 


-82.1 


0.22 








2 


-417.8 


-.%0.1 


21.6 
















19."). 


-1 :-!(). 1 


0.7 









*Aa and A(5 are the offsets from the pointing center. 
''FWHM, P.A. and rms noise are for long baselines (> 2kX). 



TABLE 3 

Pressure Normalization Values and Integrated Y{rf,oo) Values" 



Cluster Peo(SZ) Peo(X) P^o{SZ)/P^o{X) 

(10-^° ergs cm-3) 

Abell 2631 l.OOzt"-}} 1.21±" ;6 0.82±||-"a 

Abell 2204 9.90ili^;g^ 9-71111^7 1.02ifc^^ 

^sph,Sz(»'500) '5^sph,x(''500) 'i'Bph,Sz(''50o)/'5^sph,x(»'BOo) 

(10-") 

Abell 2631 9.13±;-Jg 11.13±J p 0.82±!! Jj; 

Abell 2204 44.97^^^^ 43.93ilii|^ 1.02±°;i° 



''Statistical and Chandra calibration systematics are included in the measurement of masses. 
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TABLE 4 

Results from Joint X-ray/SZ Analysis 













Model Parameters 










Cluster 


rieo 
(lO-^cm-3) 


(arcsec) 


n 







T° 
(keV) 


'"cool 

(arcsec) 


a 


7" 


Da 
(Mpc) 


Abell 2631 




26i.i±^*:« 


9.6±ii 


2.0 












799±™| 


Abell 2204 


3.90±°:?| 


22.7±ll 


6.9±l;l 


1.37±g 


10 
08 


i4.9±J:« 


2o.o±°:? 


0.17 ±0.01 


2.0 


575±tl 



Cluster Masses*^ 

Masses Evaluated At A= 2500 Masses Evaluated At A=500 



Cluster TA Mgas M^q^ /gas Afgas M^qj /gas 

(") ^^""^fa) (1"'^^q) (") (^""^ej ^^"'"^^^ 

Abell 2631 103.2±;^;^ 1.71±?:1^ 1.37±J:^^ oTiiiSo 289.3±^t6 sIsSiJ^^T 6.16±^:g 0.131±g:, 

Abell 2204 231.3±|:g 5.03±;:gf 4.09±g:|^ 0.122±g:°;^ 492.6±^g:f 13.09±i:^2 7.95±^:^^ 0.161±g:; 



^Thc reader is cautioned that To is not a global temperature, but rather a model parameter in Equation 6. 
''The parameter 7 is fixed in the model. 

"^Statistical and Chandra calibration systematics axe included in the measurements. 



TABLE 5 
Results from X-ray Analysis 













Model Parameters 










Cluster 




(arcsec) 


n 




TS 
(keV) 


rcool 
(arcsec) 


a 




Da' 
(Mpc) 


Abell 2631 


"^^0.03 


248.0±:j^-g 


9.27±; 


91 

34 


2.0 8.3±g-7 








840.0 


Abell 2204 


4 12+0.23 


22.5±;:^ 


6.76±J 


19 
82 


1.38±g:gi 15.0±gJ 


20.0±g;^ 


Q .7 ,0.08 
"■^ ' ='^0.09 


2.0 


526.0 



Cluster Masses'^ 

Masses Evaluated At A= 2500 Masses Evaluated At A=500 



Cluster 


(") 


Mgas 
(1O"M0) 


Aftot 
(10"Mq) 


/gas 


TA 
(") 


Mgas 
(IO^^^Mq) 


Mtot 

(IO^-'Mq) 


/gas 


Abell 2631 
Abell 2204 


ii3.9±^:g 

234.4±^;| 


O /f; + 0.26 
^.*OIC(, 27 


3.75±g;;« 


U.liU^O (,03 


309.9±?°-^ 


10.34±g:™ 
10.64±g:^^ 


8.oo±;-33 

7-27±g;g 


"■-^^^^0.018 
U.ITO^O (,(,7 



''The reader is cautioned that To is not a global temperature, but rather a model parameter in Equation 6. 
''Parameters 7 and Dj^ are fixed in the model. 

"^Statistical and Chandra calibration systematics are included in the measurement of masses. 



TABLE 6 
Sources of Uncertainty 







''2500 




'"500 




Effect on /gf 


IS Effect 


on Pressure 


Effect on fgs 


s Effect 


on Pressure 


Source 


(%) 


SZ (%) 


X-ray (%) 


(%) 


SZ (%) 


X-ray (%) 


Kinetic SZ effect 


±4 


±4 




±4 


±4 




Radio Point Sources 


±2 


±1 




±2 


±1 




Asphericity 


±20 


±10 


±10 


±20 


±10 


±10 


X-ray ba<;kground 


±5 




±1 


±9 




±2 


SZA calibration 


±10 


±6 




±10 


±6 




Hydrostatic Equilibrium" 


-9 






-11 






Model Assumptions 


±10 


±3 


±3 


±10 


±5 


±5 


Helium Sedimentation 


+10 




i4 


+5 






Total Systematic 


±28 


±13 


±11 


+27 
=■=29 


±14 


±12 



Chandra calibration uncertainties'* 




Surface Brightness 


±10 


Temperature 


±1 



''Uncertainty is theoretically motivated by Lau et al. (2009) as discussed in Section 4.3.4. 

''These systematic uncertainites are added to the data prior to the fit, and their effect is included in the derived masses and pressure at 
all radii. 
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APPENDIX: MCMC REPARAMETERIZATION USING A SINGULAR VALUE DECOMPOSITION METHOD 

Correlation among model parameters is a common feature of analytic models such as the beta model (Cavaliere & 
Fusco-Femiano 1976), the Vikhlinin et al. (2006) model, and the Bulbul et al. (2010) model used in this paper. The 
Monte Carlo Markov chain (MCMC) method for the analysis of X-ray and SZE data described in Bonamente et al. 
(2004) accounts for this correlation, and therefore correlation is not an issue when evaluating integrated quantities 
such as masses and Y values, and their uncertainties. Strong parameter correlation, however, may cause the MCMC 
to be inefficient in its sampling of parameter space (see, e.g., Gilks et al. 1996, page 90), requiring long chains with low 
acceptance rate because of the poor mixing. A common solution is the use of a singular value decomposition (SVD, 
e.g.. Press et al. 1992) to perform a linear transformation of the parameters to reduce the correlation among model 
parameters, and increase of the rate of acceptance in the Monte Carlo Markov chain. For the X-ray analysis of the 
Chandra data of Abell 2631 shown in Table 5, the four model parameters (rieo, Ts, n and /3) are transformed into four 
SVD parameters {svdo through svds), and the usual Metropolis-Hastings MCMC is applied to the SVD parameters. 
The accepted parameters are then transformed back to the original Bulbul et al. (2010) model parameters, for which 
we calculate integrated quantities. 

The effect of the rcparametcrization is shown in Figures 5 and 6. The strong correlation present between certain 
pairs of parameters, especially and n, is absent from the SVD parameters. The values of the correlation coefficients 




7x10 ^ 7.5x10 8x10 8.5x10 7x10 7.5x10 ^ 8x10 ^ 8.5x10 7x10 7.5x10 ^ 8x10 ^ 8.5x10 

neO neO neO 




200 300 400 500 600 200 300 +00 500 600 

Rs(orcsec) Rs(orcsec) 




Fig. 5. — Scatter plot for parameters of the Bulbul et al. (2010) model applied to the X-ray analysis of Abell 2631. Notice the strong 
correlation among some of the parameters, especially n and r^. For clarity, only every 100-th parameter in the chain is plotted. 

for the original parameters and the SVD parameters are shown in Table 7. With this reparameterization, we obtain 
an acceptance rate of approximately 30%, which is a factor of few higher than the typical acceptance rate obtained 
using the original parameters. 



TABLE 7 

Correlation Coeffk;ients for the X-ray Analysis of Abell 2631 



Parameter rieo fs n Txo 
085 -0.82 0.31 
1.00 -0.37 
-0.35 



Parameter svdo svd\ svd\ svdz 

svdo ~ -0.01 -0.14 0.21 

svdi ■■■ 0.02 -0.04 

svd2 ■ ■ ■ -0.15 

svds 



fleo 

n 
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Fig. 6.— 



Scatter plot for the SVD parameters of the same chain as in Figure 5. 



